First Results from the PHOBOS Experiment at RHIC 



Christof Roland for the PHOBOS Collaboration 

B.B.Back 1 , M.D.Baker 2 , D.S.Barton 2 , R.R.Betts 6 , R.Bindel 7 , A.Budzanowski 3 , W.Busza 4 , 
A.Carroll 2 , M.P.Decowski 4 , E.Garcia 7 , N.George 1 , K.Gulbrandsen 4 , S.Gushue 2 , C.Halliwell 6 , 
G.A.Heintzclman 2 , C.Henderson 4 , R.Holyriski 3 , D.J.Hofman 6 , B.Holzman 6 , E.Johnson 8 , J.L.Kane 4 , 
J.Katzy 4 , N. Khan 8 , W.Kucewicz 6 , P.Kulinich 4 , W.T.Lin 5 , S.Manly 8 , D.McLeod 6 , J.Michalowski 3 , 
A.C.Mignerey 7 , J.Mulmenstadt 4 , R.Nouicer 6 , A.Olszewski 2 ' 3 , R.Pak 2 , I.C.Park 8 , H.Pernegger 4 , 
C.Reed 4 , L.P.Remsberg 2 , M.Reuter 6 , C.Roland 4 , G.Roland 4 , L.Rosenberg 4 , P.Sarin 4 , P.Sawicki 3 , 
W.Skulski 8 , S.G.Steadman 4 , G.S.F.Stephans 4 , P.Steinberg 2 , M.Stodulski 3 , A.Sukhanov 2 , J.-L.Tang 5 , 
R.Teng 8 , A.Trzupek 3 , C.Vale 4 , G.J.van Nieuwenhuizen 4 , R.Verdier 4 , B.Wadsworth 4 , F.L.H.Wolfs 8 , 
B.Wosiek 3 , K.Wozniak 3 , A.H.Wuosmaa 1 , B.Wyslouch 4 

1 Physics Division, Argonne National Laboratory, Argonnc, IL 60439-4843 
2 Chemistry and C-A Departments, Brookhavcn National Laboratory, Upton, NY 11973-5000 
3 Institute of Nuclear Physics, Krakow, Poland 
4 Laboratory for Nuclear Science, Massachusetts Institute of Technology, Cambridge, MA 02139-4307 
5 Department of Physics, National Central University, Chung-Li, Taiwan 
6 Department of Physics, University of Illinois at Chicago, Chicago, IL 60607-7059 
7 Department of Chemistry, University of Maryland, College Park, MD 20742 
8 Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627 

PHOBOS is one of four experiments studying Au-Au collisions at RHIC. During the first 
running period RHIC provided Au+Au collisions at y/s NN = 56 GeV and 130 GeV. The 
data collected during this period allowed us to study the energy and centrality dependence of 
particle production, the anisotropy of the final state azimuthal distribution and particle ratios 
at mid-rapidity. 



1 Introduction 

In June 2000, the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory 
delivered the first collisions of gold nuclei at center of mass energies several times larger than 
that previously available at other accelerators. The main goal of the RHIC research program is 
to study the behavior of strongly interacting matter under conditions of extreme temperatures 
and energy densities satisfying the prerequisites predicted for the creation of the Quark Gluon 
Plasma (QGP). 

2 Experimental setup and event selection 

The PHOBOS detector employs multiple arrays of silicon pad detectors to perform track recon- 
struction, vertex detection and multiplicity measurements. The main components are: a single 
layer multiplicity detector which covers the pseudo-rapidity range \rj\ < 5.4 with nearly complete 
azimuthal coverage, a two-layer vertex detector covering \rj\ < 1.5 and a magnetic spectrometer 




Figure 1: Measured pseudorapidity den- 
sity normalized per participant pair for cen- 
tral Au+Au collisions. Systematic errors are 
shown as shaded areas. Also shown are results 
of Pb+Pb data (CERN SPS), HIJING 11 sim- 
ulations and a parametrization of pp data. 



Figure 2: Normalized pseudorapidity density 
dAT c h/d77||^| <1 /(0.5 x N par t) as a function of 
the number of participants. The solid square is 
from Ref. 1 . Theoretical calculations are shown 
from HIJING 7 (solid line), KN 8 (dotted curve) 
and EKRT 9 (dashed curve). 



with 16 layers providing momentum measurement and particle identification near mid-rapidity. 
Details of the experimental arrangement are described inlira. 

The primary event trigger was provided by two sets of 16 scintillator paddle counters which 
detected charged particles in the pseudorapidity range 3 < \r]\ < 4.5. The paddle counters were 
also used for offline event selection and determination of the collision centrality. Details of this 
procedure can be found inta 



3 Results 



Multiplicity and charged particle density. The first observable studied was the energy 
dependence of the charged particle pseudo-rapidity density near mid-rapidity for central col- 
lisions, dNch/drjU^i. The measurement is described in detail mil. We observed a primary 
charged particle density of dN c h/ dri^^i = 408 ± 12(stat) ± 30(syst) at y/s NN = 56 GeV and 
555 ± 12(stat) ± 35(syst) at y/s NN = 130 GeV. Normalizing per participant pair, we deduce 
dN ch /dr]\ M<1 /(0.5 x N part ) = 2.47 ± O.l(stat) ± 0.25(syst) and 3.24 ± O.l(stat) ± 0.25(syst), 
respectively. The data were corrected for backgrounds from decay feed-down, secondary in- 
teractions, stopping particles in the beampipe and detector inefficiencies using Monte Carlo 
simulations of our apparatus. The normalized yield per participant obtained for Au+Au colli- 
sions, proton-antiproton (pp) collisions i and central Pb+Pb collisions! are presented in Fig. [l|. 
Central Au+Au collisions show a significantly larger charged particle density per participant 
than non-single diffractive (NSD) pp collisions at comparable energies. The observed 31% in- 
crease of density from 56 to 130 GeV in central Au+Au collisions is significantly larger than the 
increase for pp data in the same energy interval!. Comparing our data to those obtained at the 
CERN SPS for Pb+Pb collisions at y / sj^" = 17.2 GeV, we find a 70% higher charged particle 
density per participant near r] = at y/s NN = 130 GeV. The dependence of the charged particle 
density dN^/drjl^^i on event centrality together with theoretical predictions Bi'i are shown in 
Fig. | 

Further information on the underlying physics can be obtained by studying the change in 
particle production over a wider range in rj. We performed this measurement using the PHOBOS 
multiplicity detector, consisting of the Octagon barrel covering <3.2 and the Ring counters 
detecting particles with 3 < ln\ < 5.4. Details of the hit-counting procedure and the background 
correction can be found in 113. For the most central 3% of the collisions we find an average 
number of charged particles of 4100 ± 100(stat) ± 400(syst) within |?7| < 5.4. 

The measured dN c h/dr] distributions scaled by 0.5 x N par t for three selected centrality bins 
are shown in Fig. ||a. For comparison, the predictions of HIJING EJ appear in Fig. ||b. At 
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Figure 3: (a). Measured dN c h/ df]/ '({N pa rt) /2) for {N par t) — 102 (circles), 216 (squares) and 354 (diamonds). We 
estimate a systematic uncertainty of 10% independent of r\. (b). Same as (a) from HIJING. 

midrapidity, the non-central data are higher than the predictions. Both data and HIJING show 
that the bulk of additional particle production in central collisions, as compared to peripheral 
ones, occurs near mid-rapidity. The central plateau observed in the data is significantly wider 
than the HIJING prediction. 

Event Anisotropy. The momentum space anisotropy, which is a reflection of the initial-state 
space anisotropy for non-central collisions, is commonly quantified by the first and second Fourier 
coefficients of the azimuthal particle distribution relative to the event-by-event reaction plane, 
dN /d{4> — ^ r) = Nq[1 + J2i cos[i(4> — ^r)]), where is the reaction plane angle. Usually, 
the first Fourier component (yi) is referred to as directed flow and the second one (i^) as elliptic 
flow. 

In PHOBOS, the elliptic flow was calculated from the azimuthal distribution of hits in the 
Octagon detector. For every event the event plane was determined by a fit to the azimuthal hit 
distribution. The observed raw value of V2 was corrected for the effects of detector granularity, 
the resolution of the event-by-event determination of the event plane relative to the true reaction 
plane and the contribution of background hits. Details of the analysis can be found inta. 

We determined V2 for events in seven classes of centrality based on the normalized paddle 
signal. We observed a strong centrality dependence of v% ranging from 0.03 for the most cen- 
tral events up to 0.07 for peripheral events (Fig. ||). We have also studied the pseudorapidity 
dependence of V2 over a large range in centrality (Fig. |5|). Outside of the central pseudorapidity 
region, i.e. \rj\ > 1, V2 drops very quickly, approximately following the shape of the charged 
particle distribution (see Fig. |3|). 

Particle ratios. The ratio of antibaryons to baryons plays a key role in understanding the 
properties of the hot and dense system formed in heavy-ion collisions at high energies. Us- 
ing the mid-rapidity spectrometer, we have studied the ratios of multiplicities of particles and 
antiparticles for primary charged pions, kaons and protons at ^s NN = 130 GeV. 

The particle ratios were measured for the 12% most central events, for which we estimate 
{Npart) = 312 ± 10(syst). As the geometrical layout of the PHOBOS detector leads to an asym- 
metry in the acceptance and detection efficiency for positively and negatively charged particles 
in the same events, we used data taken with both polarities of the PHOBOS magnet. The most 
important systematic effects such as geometrical acceptance and tracking efficiency cancel in the 
ratio. The raw particle ratios are corrected for particle decays, secondary interactions and the 
contribution from feed-down of strange hadrons. Further details can be found inlll. Within our 
acceptance we find the following ratios: 
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Figure 4: Elliptic Flow V2, in the region -1.0< r/ <1.0 
as a function of normalized paddle signal. The error 
bars are statistical only. The estimated systematic 
error is Av2 = 0.007. 



Figure 5: Elliptic Flow V2, averaged over centrality, 
as a function of pseudorapidity r\. The error bars 
are statistical only. The estimated systematic error is 
A«2 = 0.007. 



(vr-)/(7r + ) = 1.00 ± O.Ol(stat) ± 0.02(syst) 
(K-) /(K+) = 0.91 ± 0.07(stat) ± 0.06(syst) 
(p) /(p) = 0.60 ± 0.04(stat) ± 0.06(syst) 



We estimated a baryo-chemical potential [j,b from the (K~) / (K + ) and (p) / (p) ratios, using 
a statistical model calculation ta. For a realistic range of freeze-out temperatures of 160 to 
170 MeV, both ratios are consistent with hb = 45 ± 5 MeV. This value is much lower than 
Hb = 240-270 MeV obtained in thermal fits to Pb+Pb data from the CERN SPS00, showing 
a closer but not yet complete approach to a baryon-free regime at RHIC. 
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